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ABSTRACT 

We report on interferometric observations at 1.3 mm at 2" —3" resolution using the Combined Array 
for Research in Millimeter- wave Astronomy (CARMA). We identify multi- wavelength counterparts of 
three submillimeter galaxies (SMGs; Fimm > 5.5 mJy) in the COSMOS field, initially detected with 
MAMBO and AzTEC bolometers at low, - 10" - 30", resolution. All three sources - AzTEC/Cl, 
Cosbo-3 and Cosbo-8 - are identified to coincide with positions of 20 cm radio sources. Cosbo-3, 
however, is not associated with the most likely radio counterpart, closest to the MAMBO source 
position, but that further away from it. This illustrates the need for intermediate-resolution (^ 2") 
mm-observations to identify the correct counterparts of single-dish detected SMGs. All of our three 
sources become prominent only at NIR wavelengths, and their mm-to-radio flux based redshifts suggest 
that they lie at redshifts z >2. As a proof of concept, we show that photometric redshifts can be well 
determined for SMGs, and we find photometric-redshifts of 5.6 ±1.2, 1.9lo;5, ^-^d ^ 4 for AzTEC/Cl, 
Cosbo-3, and Cosbo-8, respectively. Using these we infer that these galaxies have radio-based star 
formation rates of > 1000 M0yr~^, and IR luminosities of ^ 10^^ L© consistent with properties of 
high-redshift SMGs. In summary, our sources reflect a variety of SMG properties in terms of redshift 
and clustering, consistent with the framework that SMGs are progenitors of z ~ 2 and today's passive 
galaxies. 

Subject headings: galaxies: fundamental parameters - galaxies: active, evolution - cosmology: obser- 
vations - radio continuum: galaxies 



1. INTRODUCTION 

Submillimeter galaxies (SMGs) are ultra-luminous, 
dusty star-bursting systems with extreme star formation 
rates in the range of ^ 100 — 1000 M© yr~^ (e.g. Blain 
2002). It has been shown that the bulk of this population 
is between z~2 and 3 (e.g. Chapman et al. 2005). But re- 
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cently a possible high-redshift tail of SMGs has started 
to emerge (e.g. Younger et al. 2007, 2009; Vahante et 
al. 2007). To dat e, only ab out ten z > 4 SM Gs have 
been confirmed (Da ddi et al .' 2009a, b:' Capak e t al. 200^ 
Schinnerer et al. 2008; Ricchcrs et al. 2010; Capak et aU 
2010; Smo lcic et al.. .2011.: .Coppin et al.. .2009.. .2010r 
Knudscn et al.l l2010l : iCox et all 120111: ICombes et al.l 
2012|). Their number density is still consistent with 



that expected in co s mological mod els (Baugh et al. 2005]; 
Copp in et all 12009: Sm olcic et al.l 1201 If ). Note however 
that these z > 4 studies are not complete, and may even 
point to the existence of a new or different SMG pop- 
ulation (|Wall et al.l I2008D . The intense starburst that 
creates the submillimeter bright emission is likely to oc- 
cur when the bulk of the stellar mass is being assem- 
bled in these galaxies; SMGs are generally believed to be 
the progenitors of today's massive red-and-dead elliptical 
gala xies which formed in an intense burst at high redshift 
(e.g. lCimatti et al1l2008D . It is therefore critical to study 
in detail these cosmologically important objects. 

SMGs are generally detected in mm and sub-mm sur- 
veys with single-dish telescopes that have large beams 
(> 10"). The next step is then to pinpoint the pre- 
cise locations of these objects and to match them with 
their multi-wavelength counterparts and obtain a red- 
shift. Finding the real counterpart for an SMG is not 
trivial because the spatial density of optical/IR galaxies 
in deep fields is high and usually there are multiple galax- 
ies within one single dish mm/sub-mm beam. Deep ra- 
dio, mid-IR, optical and UV (and hard X-rays for AGN) 
data of higher resolution have been used to identify the 
right counterpart by tracing the bright star formation 
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or AGN activity (e.g. Ivison et al. 2007). However, de- 
pending on galaxy properties and redshift, these differ- 
ent tracers are likely to introduce identification-biases, 
i.e. provide true identifications for only a fraction of the 
sample (and likely at low redshifts). In essence, the most 
efficient and least biased way to associate counterparts 
is through high-resolution mm observations. This has 
to date been a time-consuming process that resulted in 
a total of ^ 50 SMGs detected via mm-interferometry 
(Downes et al. 1999, Frayer et al. 2000; Dannerbauer et 
al. 2002; Downes & Solomon 2003; Genzel et al. 2003; 
Kneib et al. 2005; Greve et al. 2005; Tacconi et al. 
2006; Sheth et al. 2004, 2012, lono et al. 2006; Younger 
et al. 2007, 2009; Aravena et al. 2010; Ikarashi et al. 
2011; Tamura et al. 2010; Katsukade et al. 2010; Wang 
et al. 2011; Chen et al. 2011; Neri et al. 2003; Chap- 
man et al. 2008; Hatsukade et al. 2010). To date the 
largest comprehensive sample of SMGs detected at high- 
resolution via mm-interferometry is that in the COSMOS 
field (Younger et al. 2007, 2009, Aravena et al. 2010) and 
consists of ~ 20 sources in total. Here we present 1.3 mm 
imaging at ~ 2 — 3" resolution with the CARMA inter- 
ferometer of 3 further SMGs (Fimm > 5.5 mJy) in the 
COSMOS field originally detected in the MAMBO- and 
AzTEC-COSMOS surveys (Bertoldi et al. 2007; Aretx- 
aga et al. 2011). The counterpart association based on 
previous data was highly ambiguous due to multiple or 
faint potential radio counterparts lacking optical/NIR 
detections. We adopt Hq = 70, iljv/ — 0.3, Ha = 0.7. 

2. DATA 

2.1. COSMOS survey 

The COSMOS project is a panchromatic (X-ray 
to radio) survey of an equatorial 2n° field. The 
field ha s been observed w ith the major space- (Chan- 
dra - lElvis et al.l 120091 GALEX - [ Zamoiski eTall 



200?; HST -'Scoville et al.|[2007t |Koekemocr ct al.loOff 



Sanders ct al. 2007; 



Leauthaud et al. 2 007:; Spitzer 

libert et al.l I2OI0I: iLe Floc'h et"alll2009l: iFraver et alj 
2009() and ground- based telescopes ( Subaru, CFHT 
UKIRT. NOAO - ICamk et al.l [2007t iTaniguchi et al. 
2007'; 'Mc Cracken et al.l l2010l : VLA iSchinnerer et al ' 



2007, 201Ci, Smolcic et al., in prep.) in more than 30 
bands. Here we additionally use deep UltraVista obser- 
vations in Y, J, H, Ks bands taken between Dec/2009 
and Apr/2010 (McCracken et al., in prep.). Fractions of 
the field have been surveyed at (sub-)mm wavelengths 
(Bertoldi et al. 2007; Scott et al. 2008; Aretxaga et 
al. 2011). The three SMGs targeted here were ini- 
tially detected by MAMBO (Bertoldi et al. 2007) and 
AzTEC/ASTE (Aretxaga et al. 2011) surveys of the 
COSMOS field. The deboosted -^ 1 mm fluxes are 
13.0;iJ mJy (Aretxaga et al. 2011), 7.45 ± 1.1 mJy, and 
5.45 ±1.0 mJy (Bertoldi et al. 2007) for AzTEC/Cl, 
Cosbo-3, Cosbo-8, respectively (see also tabidet ). 

2.2. CARMA observations, data reduction and source 

detection 

We observed the three sources at 1.3 mm using the 
CARMA interferometer in a compact - D array - config- 
uration. The targets - AzTEC/Cl, Cosbo-3, and Cosbo- 
8 - were observed with 15 antennas (corresponding to 
105 baselines) in Feb. /Mar., 2009 for a total on-source 



time of 1.5, 4.0, and 5.4 hours, respectively. Weather 
conditions varied between acceptable and very good for 
1mm observations. The nearby quasar 1058-1-015 was 
observed every 15 minutes for secondary amplitude and 
phase calibration. The strong calibrator sources 3C 84, 
3C273, and 0854-1-201 were observed at least once per 
track for bandpass and flux calibration. Radio pointing 
was performed at least every 2.5 hr on nearby sources. 
The resulting total flux calibration is estimated to be 
accurate within 15-20%. 

The upper (lower) sidebands of the 1 mm receivers were 
centered at 230 (225) GHz. Each sideband was observed 
with 45 channels each 31.25 MHz wide, for a total band- 
width of 2.8 GHz (2x1406.25 MHz). For data reduction 
and analysis, the MIRIAD package was used. The flnal 
data cube obtained after flagging (and combination of 
the datasets from all runs) was collapsed along the fre- 
quency axis to obtain 1.3 mm continuum images. The 
u ~ V data were imaged with natural baseline weighting, 
leading to synthesized clean beam sizes (rms values) of 
4.2" X 3.1" (2.3 mJy/beam), 2.7" x 1.9" (0.7 mJy/beam), 
and 2.6"x2.4" (1.5 mJy/beam) for AzTEC/Cl, Cosbo-3, 
and Cosbo-8, respectively. The observations are summa- 
rized in tab:obs . 

The 1.3 mm CARMA stamps are shown in Fig.[T] All 
three sources are detected at a ^ 3 — 4ct level. We stress 
that the positions of our 1.3 mm sources perfectly coin- 
cide (within < 0.3") with significant 20 cm (1.4 GHz) 
radi o detections drawn from th e VLA-COSMOS sur- 
vey (ISchinnerer et al.ll2007l 120101 see tabidet ). As the 
chance probability of finding a radio source within the 
CARMA beam (given the radio source number density) 
is of the order of only ~ 10^**, this significantly boosts 
the validity of our mm detections. We have extracted the 
1.3 mm (230 GHz) fluxes using the AIPS tasks MAX- 
FIT (that identifies the position with maximum value in 
a selected pixel array) and JMFIT (that fits a 2D Gaus- 
sian to selected pixel arrays). The flux densities and the 
corresponding errors are summarized in tabidet . The 
fluxes are in relatively good agreement (i.e. within la) 
with those inferred from the MAMBO data, and show a 
stronger deviation from the AzTEC 1.1 mm data. This 
is likely due to a steep spectral index in the rest-frame 
sub-mm band which translates into a rapid change in flux 
even within the 1 mm window. Assuming /3 = 1.0 and 
2.0 we expect a factor between 1.6 and 1.9 discrepancy 
between the observed 1.3 mm and 1.1 mm flux densities. 

3. ANALYSIS AND RESULTS 

3.1. Multi-wavelength counterparts and their 
photometry 

Images of the three SMGs at various wavelengths, 
with 1.3 mm CARMA contours overlaid, are shown in 
Fig. [1] All three sources are coincident (within < 0.3") 
with 20 cm radio detections as summarized in tabidet . 
Only Cosbo-8 is detected in X-ray emission (within the 
Chandra-COSMOS survey; Elvis et al. 2010). The muhi- 
wavelength counterparts of our sources become promi- 
nent only at NIR wavelengths, and are blended by nearby 
sources. In order to extract the most accurate photom- 
etry for these sources we have carefully deblended the 
SMG emission (up to the Spitzer/IRAC 8 /im band), 
and extracted its flux as described in the following. 
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TABLE 1 

Summary of observations with the CARMA interferometer at 230 GHz in D-configuration. 



Source 



Pointing position 
[J2000] 



Date 



On-source 
time [hr] 



beam 

farcsccl 



rms 
[mjy/beam] 



AzTEC/Cl 10 01 41.68 
Cosbo-3 10 00 57.20 
Cosbo-8 10 00 00.00 



+02 27 11.80 
+02 20 13.00 
+02 06 34.00 



Mar/2009 
Feb/2009 
Mar/2009 



1.5 
4.0 
5.4 



4.2" X 3.1" 
2.7" X 1.9" 
2.6" X 2.4" 



2.3 
0.7 
1.5 



TABLE 2 
CARMA detections and VLA counterparts 



Source CARMA position i^l.Smm i^l.2mm Pl.lmm VLA -Fl.4GHz 

name [J2000] [mJy] [mJy] [mJy] distance [^Jy] 



-1.0 

hl.l 

-1.0 

hl.l 

-1.2 



photo-z photo-z 
mm/radio UV-MIR 



^■•^-1.4 

,+0.6 

-1.0 

1+0.5 

-0.7 



AzTEC/Cl=» 
Cosbo-3'' 
Cosbo-8'' 



10 01 41.75 
10 00 56.95 
09 59 59.92 



+02 27 13.06 
+02 20 17.79 
+02 06 33.41 



7.4 + 2.3* 
5.4 + 0.7*" 
4.8 + 1.5* 



7.45 + 1.1 
5.45 + 1.0 



13: 
9.6; 
3.7; 



0.28" 
0.25" 
0.35" 



44 + 10+ 
78 + 13+ 
104 + 13+ 



3.2: 
1.9j 



5.6 + 1.2 
1.9j 

r^ 4 



,+0.9 
-0.5 



* Peak intensity 

** Integral intensity given; the peak intensity is 2.8 mjy 

+ Adopted from the VLA-COSMOS catalogs (Schinnerer et al. 2007, 2010) 

"" Arotxaga et al. (2011) 

^ Bcrtoldi et al. (2007) 



230GHz 



3.6/^rn 



B/.im 




Fig. 1. — NIR to radio stamps for AzTEC-Cl, Cosbo-3, and Cosbo-8 with 1.3 mm (CARMA) contours overlayed. The contour levels 
are at 2a, 3cr, and 4a (la = 2.3, 0.7, 1.5 mJy/beam for AzTEC-Cl, Cosbo-3, and Cosbo-8, respectively). The cross indicates the pointing 
center at 230 GHz. The circle (2" in diameter) indicates the CARMA source position. 



3.1.1. Calibration of photometry using the COSMOS 
photometric- catalog 

In this section we describe our photometric extrac- 
tion procedure that wiU be apphed in the next section 
to the CARMA-COSMOS UV-MIR counterparts. The 
photometry of sources in the COSMOS photometric cat- 
alog is extracted using aperture techniques, which we 
also adopt here. To validate our photometric extrac- 
tion procedure, and to estimate its uncertainty we have 
drawn ~ 10 rando m sources from the i-band selected 
catalog (jCapak et al.1,2001^1 that also have IR detections, 



and are outside masked regions in the field. Using im- 
ages at their orig inal resolution (no t convolved to a com- 
mon FWHM; see lCapak et al.ll2007l for details), we adopt 
aperture sizes for individual bands as summarized in 
tab:apert . These aperture sizes were chosen to achieve 
the best agreement with the photometry reported in the 
COSMOS photometric catalog. We also corrected our 
aperture magnitudes for slight systematic offsets (see 
tab:apert ) in order to put them on the COSMOS pho- 
tometry scale. Thus, our final magnitudes are computed 
?Tiapcrt+"Zoffsct- Notc also, that for the IRAC 



as m 



final 
apcrt 
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TABLE 3 
Photometry calibration 



band 


Aperture 


mag 


mag 




diameter ["] 


offset 


error 


g 


3 


0.1123 


0.034 


r 


3 


0.2514 


0.016 


i 


3 


0.2619 


0.034 


z 


3 


0.2203" 


0.039 


J 


3 


0.0828 


0.063 


H 


2.3 


0.0644 


0.054 


K, 


2 


-0.0390 


0.080 


Y-UltraVista 


2 


0.0000 


0.040 


J-UltraVista 


2 


-0.0225 


0.045 


H-UltraVista 


2 


-0.1111 


0.035 


Ks-UltraVista 


2 


-0.1225 


0.030 


3.6 fim* 


3.8 


0.0729'' 


0.048 


4.5 /^m* 


3.8 


0.0724^= 


0.038 


5.8 Mm* 


3.8 


0.0807 


0.063 


8.0 Mm* 


3.8 


0.0823 


0.148 


IA427 


2.5 


0.0828 


0.052 


IA464 


2.5 


-0.067 


0.049 


IA484 


2.5 


0.1418 


0.067 


IA505 


2.5 


0.1105 


0.060 


IA527 


2.5 


0.1210 


0.072 


IA574 


2.5 


-0.0172 


0.042 


IA624 


2.5 


0.1786 


0.051 


IA679 


2.5 


0.0301 


0.055 


IA709 


2.5 


0.0173 


0.050 


IA738 


2.5 


0.1119 


0.055 


IA767 


2.5 


0.0263 


0.087 


IA827 


2.5 


-0.0875 


0.036 



* total (i.e. aperture corrected) magnitude 

" for the computation of the z-band magnitude we adopt the rela- 
tion presented in McCracken et al. (2010), and scale it additionally 
with the offset magnitude given here 

^ in addition, for ms g^m > 21.5 we require map = 1.215 map — 
5.356 



in addition, 
i.770 



for m3.6(jm > 21.5 we require m,a 



1.310ma 



3.6 and 4.5 fini bands we fit the faint magnitude slope 
separately (as reported in tab:apert ), and for the Subaru 
z^ ba nd we adopt the relation given bv lMcCracken et all 

We estimate the average error of our magnitudes by 
statistical propagation of the magnitude errors reported 
in the catalog and the spread of the (catalog - aper- 
ture) magnitude difference (in cases where the average 
error of the catalog magnitudes is larger than the spread 
we adopt the latter as the error of our magnitudes). In 
summary, our aperture photometry matches that in the 
COSMOS photometric catalog very well (mean offsets 
are zero) , and the average error of such extracted magni- 
tudes is estimated to be ^ 0.05 (see tab:apert for exact 
values for each band). 

3.1.2. Dehlendmg 

All three sources presented in this work are blended 
in the near/mid-IR images. We deblend the sources by 
subtracting from the corresponding image a point source 
centered at the a-priori known position of the blending 
source. The position was obtained either from the Ks 
band image (Cosbo-3 and Cosbo-8) or the pixel with 
maximum value in the 3.6 /xm image (AzTEC/Cl; see 
Fig. [21 Fig. |3] and Fig. |4]). The peak value of the point 



source to be subtracted was optimized by repeating the 
procedure multiple times until a noise-level residual re- 
mained at the position of the blending source. The 
deblending uncertainty was then obtained from the lu 
spread in the extracted magnitudes when varying the 
peak flux by ±20%. This uncertainty corresponds to 
~ 0.1 mag or better for all bands, and it is then added 
in quadrature to the photometry error estimated in the 
previous Section in order to obtain the magnitude error. 
Assuming a 2D-Gaussian point-spread function was 
satisfactory for Cosbo-3 and Cosbo-8 in all affected 
bands (see Fig. [3] and Fig. |4]), while the source blend- 
ing AzTEC/Cl appears to have internal structure in 
IRAC 3.6 and 4.5 /.tm images. Rather than modeling 
AzTEC/Cl using a more complex (and therefore more 
uncertain) model we simply isolate the blending source 
within an aperture of 2" radius, mirror it around its diag- 
onal and then subtract this from the image (see Fig. [2]). 
This will obviously cause over-subtraction in the part 
of the aperture not associated with the source of inter- 
est, while it can be assumed that the contribution of the 
blending source has been well subtracted from the source 
of interest. The blended and deblended image cutouts 
for our sources are shown in Fig. [21 Fig. [H and Fig. [H 
and the extracted photometry for each source is given in 
tab:phot . 

3.2. Photometric redshifts 

In this Section we first investigate photometric redshift 
estimates for SMGs based on a sample of eight interfer- 
ometrically observed SMGs with spectroscopic redshift 
spanning a broad redshift range of z ~ 1 — 5 present 
in the COSMOS field (Sec. [323) • Showing that photo- 
metric redshifts can be reliably calculated for SMGs, we 
then derive photometric redshifts using the same method 
for t he thr ee CARMA-COSMOS SMGs analyzed here 
(Sec. [3221). 

3.2.1. Calibration of photometric redshifts for SMGs 

Photometric redshifts are computed by fitting opti- 
mized spectral templates to the spectral energy distri- 
bution of a given galaxy, leaving redshift as a free pa- 
rameter. The redshift is then determined via a x^ min- 
imization procedure. We use the Hyper-z code to com- 
pute photometric redshifts for our SMGs with the same 
parameterization as in Wardlow et al. (2010, 2011): Ex- 
tinction assuming a Calzetti et al. (2000) law, with red- 
dening {Ay) varied from to 5, and an allowed redshift 
range of to 7. Based on ~ 30 SMGs with spectro- 
scopic redshifts, drawn from the LESS survey, Wardlow 
et al. (2011) have shown that photometric redshifts for 
SMGs derived with Hyper-z using the above mentioned 
parameterization are estimated accurately (they find a 
median offset between the spectroscopic and photomet- 
ric redshifts of 0.023 ± 0.021; see also their Fig. 1). 

In general the quality of photometric-redshifts will de- 
pend on the choice of the spectral template library to 
be fit. Therefore, in addition to the templates provided 
by Hyper-z (similar to those used by Wardlow et al.; 6T 
hereafter; see below), we also test other sets of template 
libraries. Our spectral model libraries are summarized 
as follows: 

2T: Only two - burst and constant star formation his- 
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3.6/L4m: blended 



deblended 




5 0-5 

Arc Seconds 



5 0-5 

Arc Seconds 



lor 



4.5/L4m: blended 



deblended 




5 0-5 

Arc Seconds 



5 0-5 

Arc Seconds 



5.8/L4m: blended 



deblended 



8.0/L4m: blended 



deblended 





5 0-5 

Arc Seconds 

Fig. 2.— AzTEC/Cl: IRAC 3.6 ^m (top left), 4.5 ^lm (top right), 5.8 fim (bottom left), 8.0 /im (bottom right) blended (left panel column) 
and de-blended (right column) stamps. The contour levels are at —4cr, —2cr (dotted lines), 2a, Aa (dashed lines) and 2V, i = 3,4,5,6... 
(full lines), where the rms (a) has been derived locally. The CARMA position is marked by the thick circle in the center of the stamp 
(its size matches the photometry aperture for the given band; see tabiapert ). The blending source subtracted is outlined by the circle to 
the NE. In the 3.6, and 4.5 fim bands the deblending was performed by mirroring the NE source across its SE-NW diagonal within a 2" 
aperture (which artificially induces the negative area within the aperture in the deblended images). To deblend the 5.8, and 8.0 ^m bands 
a point source was subtracted. 



tory - templates drawn from the Bruzual & Chariot 
(2003) library (and provided with Hyper-z). 

6T: Six templates provided by the Hyper-z code: burst, 
four exponentially declining star formation histo- 
ries (star formation rate oc e^*''^ where t is time, 
and T — 0.3, 1, 3 and 5 Gyr) and a constant star 
formation history. This selection of SFH/templates 
is similar to the approach used by Ilbert et al. 
(2010) to compute stellar masses with LePhare. 

M: Spectral templates developed in GRASIL (Silva et 
al. 1998; Iglesias-Paramo et al. 2007) and optimized 
for SMGs by Michalowski et al. (2010). 

For this analysis we use eight SMGs in the COSMOS 
field with available spectroscopic redshifts of counter- 
parts determined interferometrically (Capak et al. 2008, 
2010; Schinnerer et al. 2008; Riechers et al. 2010; Smolcic 
et al. 2011; Karim et al., in prep, Sheth et al., in prep). 
We compute the photometric redshifts for these SMGs as 
described above and show the total x^ distribution as a 
function of redshift in Fig.[5l The overall match between 
the most probable photometric redshift (corresponding 
to the minimum x^) and the spectroscopic redshift is 



highly satisfactory. There are no catastrophic outliers. 



For source AzTEC-3 at 



-spec 



5.3 (Capak et al. 2010; 



Riechers et al. 2010) there are two x^ minima. However, 
the low redshift peak can be disregarded given that the 
galaxy is not detected in the radio. If it were a low- 
redshift SMG one would expect a strong radio detection 
given the depth of the VLA-COSMOS survey. 

Overall, all templates reach similar solutions, and the 
best agreement (i.e. tightest x^ distribution) between the 
spectroscopic and photometric redshifts is reached when 
using the Michalowski et al. (2010) spectral templates, 
and hereafter we adopt these for our photometric-redshift 
computation. A direct comparison between the photo- 
metric (based on M templates) and spectroscopic red- 
shifts is given in Fig. HI The errors indicate the 99% confi- 
dence interval. We find a median of -0.03, and a standard 
deviation of 0.08 in the (^^phot — 2^spcc)/(l + 2spec) distribu- 
tion. We conclude that our photometric redshift compu- 
tation is accurate for SMGs as expected based on results 
from previous studies (e.g. Daddi et al. 2009; Wardlow 
et al 2010, 2011; Yun et al 2011). 



3.2.2. Redshifts for CARMA-COSMOS SMCs 
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Fig. 3. — Cosbo-3 (indicated by the thick circle in the center of the stamp) is blended by two (NE and SW) nearby sources (also indicated 
by circles). From top to bottom the rows show IRAC 3.6, 4.5, 5.8, and 8.0 ^m bands, respectively. The contour levels are at —4a-, —2a 
(dotted lines), 2cr, 4(7 (dashed lines) and 2'cr, i = 3,4,5,6... (full lines), where the rms (a) has been derived locally. The deblending was 
done by iteratively subtracting a point source (2D-Gaussian) in all four IRAC bands. 



Following the same approach as described in the pre- 
vious Section we compute the photometric redshifts for 
AzTEC/Cl, Cosbo-3 and Cosbo-8 using their extracted 
(deblended) photometry (see tab:phot ) . The results are 
presented in Fig. [7] and tabrsed . We find photometric 
redshifts of 5.6 ± 1.2, and l.OtoJ ^^ AzTEC/Cl, and 
Cosbo-3, respectively. We stress that our photometric 



redshift for Cosbo-3 is consistent (within Az ~ 0.5) with 
the recently confirmed spectroscopic redshift (Riechers et 
al., in prep), affirming both our computation of the UV- 
MIR-based photometric redshifts, and deblending tech- 
nique. 

Cosbo-8 is detected in the X-rays within the Chandra- 
COSMOS survey (Fo.5_io koV = (2.49 ± 0.65) x 



CARMA-COSMOS 



TABLE 4 
Multi-wavelength photometry 



Telescope/Band 




AB magnitude 






AzTEC/Cl 


Cosbo-3 


Cosbo-8 


Subaru/g+ 


> 26.5 


>26.5 


>26.5 


Subaru/r+ 


> 26.5 


>26.5 


26.8 ±0.3 


Subaru/i+ 


> 26.1 


>26.1 


25.8 ±0.3 


Subaru/z+ 


> 25.1 


> 25.1 


24.8 ±0.3 


UltraVista/J 


> 23.9 


24.15 ±0.19 


> 23.9 


UltraVista/H 


>22.7 


23.64 ±0.18 


23.44 ± 0.20 


UltraVista/Ks 


> 22.4 


22.80 ±0.18 


22.50 ±0.20 


Spitzer/3.6 /im 


21.12±0.11 


20.78 ±0.11 


21.20 ±0.11 


Spitzer/4.5 /im 


21.28 ±0.11 


20.48 ±0.11 


20.81 ±0.11 


Spitzer/5.8 fim 


20.83 ±0.12 


20.35 ±0.12 


20.36 ±0.12 


Spitzer/8.0 fim 


20.12 ±0.18 


19.96 ±0.18 


20.12 ±0.18 


JCMT/AzTEC-1.1 mm 


i3.b^_00g 


13.94«;1? 


14.981^3^ 


IRAM 30-m/MAMBO-1.2 mm 


- 


14.22 ±0.16 


14.56 ± 0.20 


CARMA-1.3 mm 


14.23 ±0.34 


14.57 ±0.14 


14.70 ± 0.34 


VLA-20 cm 


19.80 ±0.25 


19.17 ±0.19 


18.86 ±0.14 



Magnitudes are total magnitudes already corrected for reddening. Limits are either adopted from lUapak et aTl IJ20071 . 
extracted from the aperture flux (J, H, Ks). 



") or 



TABLE 5 
Best fit properties given by the photometric redshift computation 



Source 



template 



Hyper-z Ay 



AzTEC/Cl SMMJ030226. 17+000624.5 0.00 

Cosbo-3 LESSJ033229.4-275619 2.55 

Cosbo-8'' SMMJ131215. 27+423900.9 3.00 



+ Reddening computed by Hyper-z; note that the templates already have intrinsic reddening as defined by Michalowski et al. (2010; see 

their Tab. A3) 

" Best fit template/j4v' for both x^ minima. 



10^^^ erg/cni^/s) suggesting the presence of an AGN. 
As AGN are variable sources with featureless power- 
law spectral energy distributions, special treatment for 
photometric redshift estimates is required. In the COS- 
MOS survey Salvato et al. (2009, 2011) have optimized 
the photometric redshift computation for X-ray selected 
sources reaching an accuracy of cr^z/ii+z^ ec) ^ 0.015. 
Salvato et al. (2011) find that the Chandra-COSMOS 
source associated with Cosbo-8 (CID 838; see Civano 
et al 2012 submitted) has a photometric redshift corre- 
sponding to Zphot = 0.82 ± 0.02 best fit by a normal 
galaxy template. However, the photometry is contam- 
inated by the presence of a nearby object. Deblending 
and extracting the photometry as described here, with 
the addition of J, H, Ks from UltraVista, the photometric 
redshift value ranges from Zphot = 3.6—4.3, depending on 
the luminosity prior adopted (i.e. assuming a luminosity 
typical of a low luminosity AGN, i.e. —8 < Mb < —24, 
or typical of a QSO, i.e. —20 < Mb < —30, respectively; 
Veron & Veron 1998). Assuming a low luminosity AGN 
the best fit template is a Seyfert 1.8 drawn from the Pol- 
letta et al. (2007) library while assuming a luminosity 
typical of QSOs, the best fist template is a hybrid cre- 
ated using a ULIRG (IRAS22491) and a QSO (see Fig. [7] 
and Salvato et al 2009 for details). The solutions are 
consistent within la as in both cases the redshift proba- 



bility distribution function shows a broad range of pos- 
sible solutions, rather than a well defined unique peak. 
These solutions are also consistent with the photometric 
redshift value computed as described in the previous sec- 
tion which yields two x^ minima (at z ^ 3 and z ~ 4; see 
Fig. [7]). Comparing the results from the various spectral 
libraries the best x^ value is obtained when using a QSO 
prior yielding a redshift of Zphot = 4.11q;5 (where the 
error is a la error; see Fig. [7|). Hereafter we take z ^ A 
as the best redshift estimate for this source, noting that 
spectroscopic redshift follow-up is required to disentan- 
gle between the various photometric redshift solutions 
for this SMG. 

For comparison, the mm-to-radio fiux ratios of the 
sources, regularly utilized as a redshift estimate for SMGs 
(Carilli & Yun 2002), suggest that all three sources lie at 
high (z > 2) redshift (see tab:det ). An improved version 
of the dust-independent Carilli & Yun (2002) redshift 
estimator via the observed mm-to-radio flux density ra- 
tio (Yun et al. 2011) yields redshifts of 4:.5tll, ^-"^-i.o^ 
and 1.9l[5;^ ^r AzTEC/Cl, Cosbo-3, and Cosbo-8, re- 
spectively (using AzTEC 1.1 mm measurements for all 
sources). 

In summary, we find photometric redshifts of 5.6 ± 1.2, 
1.9;°J, and - 4 for AzTEC/Cl, Cosbo-3, and Cosbo-8, 
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Fig. 4. — Blended and deblended Cosbo-8 images. H, Ks, and IRAC bands (indicated in each panel) needed to be deblended. The 
contour levels are at — 4(t, —2cr (dotted lines), 2a, icr (dashed lines) and 2'(7, i = 3,4,5,6... (full lines), where the rms {a) has been derived 
locally. The CARMA position is marked by the thick circle in the center of the stamp (its size matches the photometry aperture). The 
source that was subtracted (assuming a 2D-Gaussian point-source model) is outlined by the circle to the NE. 
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Fig. 5. — Full x^ distribution of the photometric redshift determination for SMGs in the COSMOS field with interferometrically determined 
counterparts with spectroscopic redshifts. The photometric redshifts were derived using various sets of spectral models (see text for details): 
2T (dotted lines), 6T (dashed-lines), M (full lines). The spectroscopic redshifts are indicated by vertical lines. Note that for AzTEC-3 the 
radio non-detection rules out the first x^ minimum. 

CARMA detection{3 The deblended images and pho- 
tometry for AzTEC/Cl are presented in Fig. [5] and 
tab:phot . 

The mm-to-radio flux ratio of AzTEC/Cl yields a red- 
shift of ^mm/radio = 'i-'i-i 4: Consistent with the photo- 
metric redshift of 2;phot = 5.6 ± 1.2 derived from its UV- 
MIR photometry. The inferred high redshift is consistent 
with the source not being detected at wavelengths shorter 
than 3 /im. At such a redshift (4.3-5.6) the radio flux den- 
sity of the galaxy (F20cm = 44 ± 10 /xJy) would imply a 
20 cm luminosity of (6 - 10) x lO^^ W Hz"!. If the entire 
radio emission arises from star formation in the galaxy, 
and if at these redshifts locally determined radio-star for- 
mation rate calibrators (Bell et al. 2003, Yun et al. 2001) 
can be applied (as would be suggested by the constancy 
of the FIR-radio correlation out to high redshifts; Sargent 
et al. 2010a, 2010b, Murphy 2009), this radio luminosity 
would imply a SFR of ~ 3200 - 5600 Mq yr'^ This 
is somewhat in excess of expectations for typical SMGs, 
thus it may be possible that part of the (radio) emis- 
sion from this source arises from black hole accretion. 
However, it is worth noting that such properties are not 
unusual for z > 4 SMGs. For example, the properties of 
AzTEC/Cl are very similar to those of AzTEC-1 - the 
brightest SMG in the AzTEC/JCMT COSMOS survey 
(Scott et al. 2008; Younger et al. 2007, 2009, Smolcic et 
al. 2011;F20cm = 42^Jy, z = 4.6). 




Fig. 6. — Comparison of spectroscopic and photometric redshifts 
for eight SMGs from the COSMOS field using the Michalowski et 
al. (2010) spectral templates. The errors show the 99% confidence 
interval. The median offset and standard deviation of the Az/{1 + 
^spcc) distribution are indicated in the bottom panel. 



respectively. Below we summarize the properties of each 
SMG. 

3.3. Properties of individual sources 

3.3.1. AzTEC/Cl 

AzTEC/Cl has neither a J, H or Ks band counter- 
part , while it can be associated with a source at 3.6 /Ltm 
(and becoming most prominent at 8.0 y^m) that is, how- 
ever, strongly blended with a bright source 1.97" to 
the NE (see Fig. [1] and Fig. [2]). The radio-counterpart 
of AzTEC/Cl is only 0.28" away from the reported 



^^ Based on the S/N = 4.4 radio detection at a resolution of 
O = 1.5" the expected astrometric accuracy is 0/(S/N) = 0.34" 
(note that the overall VLA-COSMOS astrometric accuracy is esti- 
mated to be better than 0.130"; see Schinnerer et al. 2007 for de- 
tails). In order to assess the astrometric accuracy of our CARMA 
detection we imaged 3C273, our secondary calibrator ~ 37deg 
away from the AzTEC/Cl field phase center. We recover its po- 
sition within 0.11" of the nominal position. This yields that the 
positional uncertainty of AzTEC/Cl is likely better than this value, 
and thus it rules out the possibility that the bright IR galaxy ~ 2" 
away from AzTEC/Cl (which also corresponds to the closest source 
to AzTEC/Cl detected in the optical) is its counterpart (or that 
of the radio source) . 
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Fig. 7. — Top panel: Total x^ distribution as a function of redshift for our SMGs (indicated in the panels). Solid curves show the 
solution based on the Michalowski et al. (2010) spectral model library used to derive the photometric redshifts as described in Sec. 13.21 For 
Cosbo-8 (right panel) we also show the total x^ distribution when using model libraries with AGN templates and a low-luminosity AGN 
prior (dashed line) and a QSO prior (dotted line) are used (see text for details). Bottom panel: The spectral energy distribution for our 
SMGs (indicated in the panels; dots and arrows for 5cr upper limits) with the best fit spectral model (corresponding to that yielding the 
minimal x^ value) overplotted (black, gray and rod lines show the best fit Michalowski et al., AGN and QSO library spectra, respectively). 
The best fit redshift is also indicated. 



3.3.2. Cosbo-3 

A source coincident with the position of Cosbo-3 is 
detected in J-band, as weh as in longer wavelength bands. 
In the Spitzer images two surrounding sources (1.9" to 
the NE and 2.4" to the SW, respectively) are blending 
its IR emission (see Fig. [1] and Fig. [3]). The deblended 
images and photometry for AzTEC/Cl are presented in 
Fig. [3] and tab:phot . 

Within the MAMBO 11" beam there are 2 radio 
sources present, at separations of 1.3" and 5.9", re- 
spectively. Contrary to expectations, the Cosbo-3 mni- 
source identified by CARMA is coincident with the NW- 
radio source and not the radio source (at Zphoto = 
2.4) closest to the mm-source identified by MAMBO 
(jBertoldi et al.ll2007t ). Although consistent (within la) 
with the MAMBO 1.1 mm fiux, the CARMA 1.3 mm 
flux density is somewhat lower. Thus, it may be possible 



that Cosbo-3 at ~ 11" resolution is itself a blend of two 
mm-sources, one of which was not detected within the 
CARMA 1.3 mm sensitivity. Our CARMA observations 
put a 3cr upper limit to the emission of a potential second 
mm-source of 2.1 mJy at 1.3 mm. 

Cosbo-3 was found to be located in a strong overden- 
sity (30 X higher than the field) of star forming galaxies 
(Aravena et al. 2010). Thus, it is possible that part of 
the MAMBO emission is distributed over several sources, 
consistent with our CARMA observations. All the galax- 
ies in the overdensity have photometric redshifts in the 
range z = 2.2 — 2.4, providing strong statistical support 
to the photometric redshift of our identified counterpart. 

The mm-to-radio flux ratio suggests a redshift of z = 
3.2 



+0.6 



Our photometric redshift, based on the de- 
blended UV-MIR data, yields Zphot = l-Qto.^g, which 
is consistent with the source's spectroscopic redshift 
(Riechers et al., in prep.) and closer to the photometric 
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redshift of the surrounding overdensity. Assuming the 
UV-MIR based photometric-redshift value the source's 
20 cm radio flux density (i^20cm = 78 ± 13 /iJy) implies 
a radio luminosity of '^ 2 x 10^** W Hz~^ and a SFR of 
~ 900 Moyr~^. Scaling an Arp 220 template we find 
an IR luminosity of Lir ~ 1.5 x 10^^ Lq. The redshift 
and star formation rate of Cosbo-3 are fairly typical for 
SMGs, found to form stars at similar rates and populat- 
ing the redshift range z = 2 — 3 (e.g. Chapman et al. 
2005; Wardlow et al. 2010; Yun et al. 2011). 

3.3.3. Cosbo-8 

A source coincident with the position of Cosbo-8 is 
detected in the radio band at high significance (i^20cm = 
104 ± 13 /iJy). As in the case of the other two SMGs, 
in Spitzer images it is blended with a source 1.2" to the 
NE (see Fig. [T] and Fig. |3]). Its deblended photometry 
is presented in tab:phot and the deblended images in 

Fig. a 

The mm-to-radio flux ratio suggests a redshift of 2; = 



1.9 



+0.5 

-0.7- 



The UV-MIR-based photometric redshift is in 
the range of 2; ^ 3.6 — 4.3 (using spectral models typi- 
cal for AGN as this source is detected in the X-rays by 
Chandra). Assuming z ^ 4 the radio flux of Cosbo-8 
implies a 20 cm luminosity of ^ 10^^ W Hz^^ and a 
radio-based SFR of ^ 6400 Moyr~^, an IR luminosity 
of LiR ^ 1.1 X 10^^ L© (based on a scaled Arp 220 tem- 
plate), and an IR-based SFR of - 1700 Moyr"! (BeU 
2003). The difference in the radio- and IR-based SFRs 
suggests the presence of an AGN also at radio wave- 
lengths. 

Cosbo-8 corresponds to a point source detected at high 
significance in both the full C-COSMOS and the best- 
PSF C-COSMOS datasets with 15 (2.3 estimated back- 
ground) and 5 (0 background) counts in the 0.5-2 keV 
band, respectively. Using the countrate-to-flux conver- 
sion factors from Puccetti et al. (2009), we obtain a flux 
of the source in the 0.5-2 keV band of (3.8 ± 1.1) x lO'^^ 
ergs s~^ cm^. Assuming z ^ A and F = 1.4, this corre- 
sponds to a rest-frame luminosity in the 2-10 keV band 
of 6.8 X lO''^ ergs s~^, which corresponds to the level of 
the emission of a typical AGN (using F = 2 results in a 
10% downward revision of the luminosity). 

In the MAMBO-COSMOS area to-date only one other 
SMG (Cosbo-11) has been confirmed as a X-ray detected 
AGN (type-1 QSO at z^poc = 1.83; Aravena et al. 2008). 
Cosbo-11 is likely on-going a merger and shows radio and 
IR luminosities consistent with purely SF activity. Based 
on the duality of properties showing both properties of 
QSO and starbust it has been classified as a starburst-to- 
QSO 'transition' system. Cosbo-8, on the other hand, is 
also on-going a major starburst as implied by its IR lumi- 
nosity, however, it shows an excess of radio emission with 
respect to that expected from the IR SED. This suggests 
that the AGN in this case is having a more important 
role in the bolometric output. Thus this source could 
also be classified as a starburst-QSO transition object, 
but possibly in a more advanced stage when the AGN 
starts to dominate the SED. 

4. SUMMARY AND DISCUSSION 

In order to unambiguously determine the multi- 
wavelength counterparts of three Fimm > 5.5 mJy SMGs 



in the COSMOS field (initially detected with MAMBO 
and AzTEC bolometers at low, > 10", resolution), 
we performed interferometric observations at 1.3 mm 
and ~ 2" - 3" resolution using CARMA. The obser- 
vations yielded 3 — 4cr detections coincident with po- 
sitions of 20 cm radio sour ces (VLA-COSMOS survey; 
iSchinnerer et al.ll2007ll2010[) . Although ah three sources 
are coincident with radio detections, our observations il- 
lustrate the need for high-resolution mm-imaging to de- 
termine the correct counterparts of bolometer-identified 
SMGs. Without high-resolution mm observations, the 
counterpart of Cosbo-3 would have been misclassified 
as our observations associate this SMG with the radio 
source (out of two radio-sources) within the MAMBO 
beam that is further away from the MAMBO source cen- 
ter. 

All three sources identified here are blended in the MIR 
by nearby bright galaxies. We have carefully deblended 
their photometry, and derived photometric redshifts. We 
find photometric redshift of Zphot(AzTEC/Cl) = 5.6 ± 
1.2, 2;phot(Cosbo - 3) == 1.9to.5, ^nd Zphot(Cosbo - 8) -- 
4. These are consistent with mm-to-radio-fiux based esti- 
mates for AzTEC/Cl and Cosbo-3 (4.3J:?:^ and 3.2t°:[;, 
respectively), but inconsistent with that inferred for 
Cosbo-8 (1.9^q'7). This is naturally understood as (part 
of) the radio fiux in Cosbo-8 may arise from the asso- 
ciated AGN identified via the X-ray Chandra detection. 
An increased radio flux due to processes not related to 
star formation would lead to an artificial decrease in the 
mm-to-radio-flux based redshift. 

Our three SMGs seem to show a relatively large red- 
shift spread, comparable to optically selected SMGs, but 
with a potential bias towards higher redshifts as all 3 
SMGs have been found to be at z > 2 (e.g. Chap- 
man et al. 2005). In general, although it has been 
shown that the SMG population peaks between red- 
shifts 2 and 3 (e.g. Chapman et al. 2005; Wardlow et 
al. 2011; Yun et al. 2011), their exact redshift distri- 
bution (and thus their cosmic evolution) is still rather 
poorly understood. This is mainly related to statis- 
tical counterpart selection biases induced by the large 
single-dish mm-beams (see e.g. Yun et al. 2011 for a 
more detailed discussion). This can be avoided by 
mm-interferometric imaging at intermediate/high angu- 
lar resolution of complete samples of SMGs. However, 
generating such samples has been a very time-consuming, 
and largely unfeasible process, and assembling complete 
samples of SMG counterparts (and their redshifts) will 
require surveys with facilities such as ALMA and LMT. 
Nonetheless, existing mm-interferometric observations of 
SMGs already suggest that a fraction of these sources 
(at least at the bright end) is unexpectedly at red- 
shifts z > 4 (e.g. 4/17 AzTEC/JCMT SMGs detected 
by SMA/VLA are spcctroscopically confirmed to be at 
z > 4; Scott et al. 2008; Younger et al 2007, 2009; Ca- 
pak et al. 2008, Schinnerer et al. 2008; Riechers et al. 
2010; Capak et al. 2010; Karim et al., in prep). In this 
cosmic epoc h to-date th e re are onl y about ten SMGs 
confirmed (Daddi et al. 2009a, b; Caoak et al. 20081 
Schinn erer et al. 200 8 ; Riec h ers ct al. 2010 ; Capa k et al ] 
2OIOI: ISmolcic et all [20TI ICoppin et al.| [20091 [201( 



■ ■- .... - 3P1" , _ 

Knudsen et all 120101: [Cox et all 120111: iCombes et al.l 
2012f ). Given their large star formation rates, these very 
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high redshift SMGs are considered to be optimal can- 
didates for the progenitors oi z ^ 2 massive red galax- 
ies (Cimatti et al. 2008). However, as this population 
is just starting to emerge their role in galaxy evolu- 
tion is still largely unexplored, and efforts to identify 
such sources and characterize their properties is critical. 
Based on min-/radio-interferometry we have associated 
AzTEC/Cl, the brightest SMG in the AzTEC/ASTE 
COSMOS survey, with a MIR/radio source at z > 4. 
Interestingly, its high- redshift, mm- and radio- fluxes are 
comparable to that of AzTEC-1 - the brightest SMG in 
the AzTEC/JCMT COSMOS survey (Scott et al. 2008) - 
with properties resembling those expected for the progen- 
itors of compact massive red galaxies at z ~ 2 (Smolcic 
et al. 2011). 

Studies of SMGs, and the formation of passive red 
galaxies (often found in the most massive galaxy clus- 
ters) suggest that the first may be progenitors of the 
second (e.g. Michalowski et al. 2010; Hickox et al. 2011). 
Consistent with this picture, in which strong clustering 
of SMGs is expected, spatial clustering analysis of z < 3 
SMGs find that SMGs cluster strongly (e.g. Hickox et 
al. 2011), and that they can be statistically associated 
with galaxy overdensities (Aravena et al. 2010). Further- 
more, only recently have two z > 4 protoclusters host- 
ing SMGs been identified (Daddi et al. 2009; Capak et 
al. 2010), providing valuable laboratories to study struc- 
ture formation at the earliest cosmic times. Based on 
i?zi^-selected galaxies Aravena et al. (2010) identified 
a significant galaxy overdensity at z ^ 2.3 in the area 
around Cosbo-3. The mm/radio- interferometric obser- 
vations and photometric redshift computation, presented 
here, strengthen the case that Cosbo-3 is indeed associ- 
ated with this galaxy overdensity (rather than being a 
fore-/back-ground galaxy). This system therefore proves 
interesting for further studies of the dense environment 



of SMGs at the peak epoch of this population, linking 
the properties of z > 4 proto-clusters hosting SMGs and 
local galaxy clusters. 

In summary, we have identified the counterparts of 
three single-dish detected SMGs in the COSMOS field 
via mm-interferometry, and presented their (deblended) 
UV-MIR photometry, and redshift estimates. Such stud- 
ies are an important step towards reaching unbiased sta- 
tistical samples of SMGs with accurately determined 
counterparts and redshifts - a necessary but yet un- 
accomplished prerequisite for comprehensive studies of 
physical properties of SMG population(s) and their role 
in galaxy formation and evolution. 
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